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Abstract 

A new dynamical solution for a gas sphere under self-gravity is 
presented to describe a development of a gas sphere from a motion-less 
state to a state of expansion with a constant speed and a reflection phe- 
nomenon in the dynamics of the surface of the sphere. 
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1 Introduction 

The polytropic model of stellar structure is often introduced in studies of 
the different stages of stellar evolution. The equation governing the static equi- 
librium of a polytropic gas sphere is called the Endem equation, whose solution 
representing the structure of star is known as the Endem solution pQ. Quasi- 
statical change of stellar structure expressed by the Endem solution is believed 
to describe the evolution of star when the Endem solution is dynamically stable. 
The dynamical stability is not guaranteed unless the polytropic index N is less 
than 3 A gas sphere of N = 3 is borderline to the dynamical stability and 
is also interesting since the polytropic index 3 corresponding to the adiabatic 
constant 7 = 4/3. There are some interesting stellar gases with 7 = 4/3, such 
as a radiation pressure dominant gas, an extremely relativistic gas and a gas 
with degenerate electrons 0- 

In order to extend the static Endem solution to dynamical states, Munier 
and Feix derived an explicit self-similar solution which describes a decelerating 
expansion of a gas sphere of N = 3 [Sj- Recently, another explicit self-similar 
solution has been obtained to describe a non accelerating nor decelerating ex- 
pansion (or contraction) of a gas sphere of N = 3 3 . 

In this paper, we present a new dynamical solution for a gas sphere of N = 3, 
which is not self similar in general but includes the three known Endem-type 
solutions described above as special cases. The new dynamical solution describes 
such interesting phenomena as a development of a gas sphere from a motion-less 
state to a state of expansion with a constant speed and a reflection phenomenon 
in the dynamics of the surface of a sphere. 
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2 Dynamical solution of gas sphere 



The equations governing the spherically symmetrical flow of a polytropic gas 
of adiabatic index 7 under the influence of its own gravitation are 

p t + up r + p(u r + ^u^j = 0, (1) 

Ut + uu r + p 1 ^ 2 p r + o> = 0, (2) 
2 

a rr H — ay = p, (3) 
r 

where the pressure P, density p, radial velocity u, gravitational potential a and 
radius r are normalized by Pq, po, uo = (^fPo/po) 1 ^ 2 , lPo/po, r o, respectively; G 
is the gravitational constant and P — p 1 . The acoustic time scale t^/uq , where 
?'o = 7-Fb/(47rGpo)j is set to equal to the time scale of free fall (4irGpo)~ 1 / 2 . 
The system of Eqs. Q, and © admits a following Lie point symmetry 



V = (t + i )fit - (7 - 2)rd r - 2pd p - (7 - 1)11$, - 7 <^, (4) 

where = oy and to is an arbitrary constant. The symmetry Q is derived 
through the standard procedure. Solving the Lie equation associated with the 
infinitesimal generator (J3J, we introduce a new independent variable T defined 
as 



dT = (5) 



and scaled variables 



_ R(x,T) U(x,T) A _ <f>(x,T) 



where 

x = -^, r = io g (t + to). 

Let us make the following ansatz 



(7) 



where 



f2(ar, T) = R(y) exp | ^ ^2 - 3[/j dT j , 
[/(x,T)-f7(T)x, (8) 
$(ar, T) - %) exp j | f 7 - 2f7J <ff j , (9) 



// = /rxp <| / ( 2- 7 - Z/Jcff 
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Figure 1: z vs. the acceleration parameter a 



This ansatz ensures that Eq. (JTJ is satisfied automatically. For 7 = 4/3 , Eqs. 
(0) and JSJl lead to 

ay + Pp~ 2 R y + $ = 0, (10) 



I7 T +I7(I7-1) = aexp<j / I 2-3U]dT (11) 



= (12) 

where a is an arbitrary constant. Eqs. i|10[) an d l|12() give a modified Endcm 
equation for the density profile R(y) 

0y> y i + ^8 y ,+8 3 + 3a = 0, 6 = R 1 / 3 , y' = -^=. (13) 

y V3 

According to numerical integration of Ea. p3|l with the boundary condition 
9(0) = 1, 6y'(0) = 0, the density profile has the first zero at finite y' = z 
for clq > —0.00219 as shown in Figure. 1. The value of z determines the radius 
of a gas sphere. For a = 0, z takes the Endem's static value z — 6.89685. Since 
the term ay' comes from the inertia term and represents acceleration effect for 
a > 0, we call a as the acceleration parameter in this paper. In the acceleration 
case, the surface of a gas sphere feels inward inertia force and shrinks from the 
Endem's static value. For a < 0, the surface feels outward inertia force and 
expands. At a w —0.00219, z diverges to infinity and a localized gas sphere 
disappears. 

It is easy to see that Ea. (|ll|l takes the following integrable form 

fu = p (14) 

where 

/(t)=exp(/ Udr). (15) 



Eq. (fH|) can be integrated once and yields 



\ft + j = b, (16) 

where b is an integration constant. In terms of f(t), the solution JHJ) is rewritten 
as 

r R{y) f , Hv) n7 , 

y = j> p = -j^ri> » = /*». ^ = t^' 

where -R(y) is expressed by a localized solution of Eq. l|13fl while Ea. (|10|l gives 
$(y) in terms of R(y). It is confirmed that the solution (|17f) satisfies the bound- 
ary condition at the moving surface of sphere. While the surface is moving with 
a speed u = \Z3f t z, the flow speed at the surface is given by the same one, i.e. 
u = fty = Viftz. 

When b = and a < 0, Eq.^TJI yields 



/w = (^( t - t .)) 2/a 

and the solution (|17|) becomes identical with the known self-similar solution |2j . 
In the absence of the inertia term ( i.e. a = 0), Ea. I|13fl reduces to the Endem 
equation. For 6 = and a — 0, f(t) is constant and Eq. l|17|l gives the static 
Endem solution. For b > and a = 0, Eq.fllBJ gives f(t) = Vbt and Ea. l(T7|) 
reproduces the expanding Endem solution which is expanding with the constant 
acoustic speed jSJ. For the other values of a and b, a new class of expanding 
solution with positive or negative acceleration is obtained. Since u = ^/3f t z, 
the acceleration rate of the surface is given by u t = a\f2>zj f 2 , of which sign is 
the same as the o's sign. 

To demonstrate the evolution of the solution l|17fl , we consider an initial 
value problem for Ea . (|16|l such that f(t)\t=o = 1 in order that y\t=o = r. 

First of all, let us consider the case in which a > and the acceleration 
rate of the surface is radially positive. Since / ^ oo as t ^ oo, ft —> V2b 
asymptotically and the surface eventually expands with a constant speed \/6bz. 
In this case, we can take a motionless initial condition so that ft\t=a = or 
u\ t= o = by choosing the arbitrary (positive) constant b as b = a(> 0). The 
expansion speed of the surface (u — y/3zf t ) is depicted for a = b = 0.01 in Fig. 2. 

Next, we consider the radially decelerating case (i.e. a < ), where we take 
ft\t=o > as an initial condition so that the initial velocity of the surface is 
positive. Then, the initially expanding surface gradually decelerates due to the 
negative inertia term. In this case, Ea. (|16fl reads 



\f? = *+^- (18) 

For b > 0, there are no reflection points and the expanding speed of the surface 
decreases to a constant speed , that is, 
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Figure 2: Expansion velocity of the surface for a=b=0.01. 
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Figure 3: Velocity of the surface for a=-0.002 and b=0.001. 



lim u = V&bz, lim u t = 0. (19) 

t — >oo t — »oo 

We illustrate the velocity of the surface for a = —0.002 and b — 0.001 in Fig. 3. 
When b < 0, a reflection point appears at / = a/b. The deceleration of the 
surface stops at the reflection point and then the gas sphere began to shrinks. 
The reflection time t r is given as When a < and b < 0, the equation l)16|) 
yields 



y/(2b-2a) a 



2b 



: arctan 



2b- 2a 



-2b 



In this reflection case, the velocity of the surface is depicted for a = —0.002 
and b = -0.001 in Fig.4, where t r = 57.485. 



5 




Figure 4: Velocity of the surface for a=-0.002 and b=-0.001. 

3 summary 

We present a new dynamical solution for a gas sphere under self-gravity which 
not only unifies the three Endem-type solutions for the polytropic index N = 3 
but also describes the following interesting phenomena. 

For positive values of acceleration parameter a, the radius of the gas sphere 
decreases from the Endem's radius due to inward inertia force, while the gas 
sphere expands for the negative a and the distinct surface disappears when 
a < —0.00219. The new solution also describes acceleration of the surface of 
a gas sphere from a motion-less state to a state of expansion with a constant 
speed for the positive a and a reflection phenomenon in the dynamics of the 
surface of a sphere for the negative a. 

References 

[1] S. Chandrasekhar, An Introduction to the Study of Stellar Structure, Dover, 
(1957). 

[2] For example, C. Hayashi, and S. Hayakawa, Astrophysics, Iwanami, (1973). 
[3] A. Munier, and M. R. Feix, ApJ 267(1983) 344. 

[4] P. J. Olver, Applications of Lie Groups to Differential Equations, Springer, 
(1986). 

[5] S. Murata and K. Nozaki, JPSJ 71 (2002) 2825. 



6 



